Second- and First-Order Phase Transitions in CDT 



J. Ambj0rn a , S. Jordan b , J. Jurkiewicz c and R. Loll 



a The Niels Bohr Institute, Copenhagen University 
Blegdamsvej 17, DK-2100 Copenhagen 0, Denmark, 
email: ambjorn@nbi.dk 

b Institute for Theoretical Physics, Utrecht University, 
Leuvenlaan 4, NL-3584 CE Utrecht, The Netherlands, 
email: s.jordan@uu.nl, r.loll@uu.nl 

c Institute of Physics, Jagellonian University, 
Reymonta 4, PL 30-059 Krakow, Poland, 
email : j urkic wicz @ th. if . uj . edu. pi 

d Perimeter Institute for Theoretical Physics, 
31 Caroline St. N., Waterloo, Ontario, Canada N2L 2Y5. 
email: rloll@perimeterinstitute.ca 



b.d 



Abstract 

Causal Dynamical Triangulations (CDT) is a proposal for a theory of quantum 
gravity, which implements a path-integral quantization of gravity as the contin- 
uum limit of a sum over piecewise flat spacetime geometries. We use Monte 
Carlo simulations to analyse the phase transition lines bordering the physically 
interesting de Sitter phase of the four- dimensional CDT model. Using a range of 
numerical criteria, we present strong evidence that the so-called A-C transition 
is first order, while the B-C transition is second order. The presence of a second- 
order transition may be related to an ultraviolet fixed point of quantum gravity 
and thus provide the key to probing physics at and possibly beyond the Planck 
scale. 
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1 Introduction 



The relation between general relativity and quantum physics is still far from being 
understood, despite decades of research. Expanding the gravitational field around 
a fixed background geometry results in a perturbatively nonrenormalizable theory. 
The best one can do within this conventional setting is to view the quantum field 
theory as an effective field theory up to a certain energy scale [TJ [2]. However, 
because of their extreme smallness, it remains unclear whether any quantum 
corrections one can compute in the range where this effective framework is reliable 
can be related to observable quantities, even in principle [3]. 

In order to define a theory which is UV-complete one has to go beyond or- 
dinary perturbative quantum field theory, either by modifying or extending the 
notion of "quantum fields" (as happens in supergravity, string theory or non- 
commutative field theory, say), or by trying to define a quantum field theory of 
gravity nonperturbatively (as, for example, in lattice or loop quantum gravity). 
In this article we will be following a nonperturbative route to quantization via the 
method of Causal Dynamical Triangulations (CDT), which falls into the latter 
category. 

One explanation for why quantum gravity could be well-defined nonperturba- 
tively, despite its lack of perturbative renormalizability, is given by the so-called 
asymptotic safety scenario. Its general framework was originally formulated by 
Weinberg jl], and is based on the hypothesis that the quantum field theory of 
four-dimensional gravity possesses a non-Gaussian ultraviolet fixed point under 
the flow of the renormalization group. The existence of such a UV fixed point 
is suggested by a (2 + e)-expansion in the dimensionality of spacetime, starting 
from the observation that two-dimensional "gravity" is formally a renormalizable 
theory [5], an observation later corroborated in [B] (for a summary of these results, 
see [7]). Further support, which does not rely on perturbing around d = 2, comes 
from using functional renormalization group methods, first applied in the context 
of gravity in a seminal paper by Reuter [8]. Explicit numerical computations of 
the renormalization group trajectories require a truncation of the infinite tower of 
higher-order curvature terms in the gravitational Lagrangian. By now, numerous 
truncations have been studied, which all point to the existence of a fixed point 
[9] (see [10] for a recent review). 

Independent evidence for the asymptotic safety of gravity - if it is realized 
- could come from a nonperturbative path integral quantization of gravity. In 
order to make the path integral well defined one usually discretizes the space of 
field configurations, which for pure gravity is the space of all spacetime geome- 
tries. This can be conveniently achieved by representing spacetime geometries 
in terms of piecewise flat, triangulated manifolds with distributional curvature 
assignments, leading to the simplicial approaches to quantum gravity, such as 
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quantum Regge calculus and dynamical triangulations (DT) (see [TT] for a re- 
view). In what follows we will study an improved, Lorentzian version of the 
latter, known as Causal Dynamical Triangulations (CDT) (see [12] for reviews). 

A crucial and nontrivial requirement for all nonperturbative approaches to 
quantum gravity is that they must reproduce well-established classical behaviour 
of gravitational physics in the low-energy limit. In this respect CDT has had some 
remarkable successes. Numerical simulations of the model have demonstrated the 
dynamical emergence of a four-dimensional classical universe as the ground state 
of the system, even though no background was ever put in by hand [14]. In 
addition, a more detailed analysis revealed that the effective action for the scale 
factor of the universe can be matched to a minisuperspace action describing a 
quantum de Sitter spacetime, including linearized quantum fluctuations [T5] . By 
contrast, the present piece of work will focus on the short-distance behaviour of 
quantum gravity, defined through CDT. 

In the following section we will recap some essential ingredients of the CDT 
approach, to set the stage for the subsequent discussion of its phase structure. In 
Sec. 13] we recall some properties of the phase diagram and describe our method 
of locating the phase transitions quantitatively. In Sec. H] we turn to the deter- 
mination of the order of two of the phase transitions, the transition line between 
phases A and C, and between phases B and C. The fact that the B-C transition 
appears to be of second order has been announced recently in [16]. Finally, our 
conclusions are presented in Sec. El 

2 Causal Dynamical Triangulations 

We begin by summarizing some important aspects of CDT, which are relevant 
for our subsequent discussion of its phase structure, and recommend consulting 
[T7] IT8] IT9] for more detailed information. The primary idea behind CDT is to 
quantize gravity in the path-integral formalism, and is best illustrated in terms 
of the analogous quantization of the nonrelativistic particle. In order to properly 
define the path integral for the latter one can discretize time by introducing a unit 
time interval of length a, and consider only paths which consist of a contiguous 
sequence of linear segments, one for each time interval. Continuum quantum 
physics is then recovered by taking the limit as a — > of the regularized path 
integral over the ensemble of piecewise linear paths. 

When generalizing from nonrelativistic particle physics to quantum gravity, 
the sum over piecewise linear (read: piecewise straight) particle paths becomes 
a sum over piecewise linear (read: piecewise flat) spacetime geometries. At this 
point one needs a concise definition of the (regularized) space of geometries to 
be summed over. The distinguishing feature of CDT, compared to other, similar 
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approaches, is that all geometries in its configuration space are equipped with a 
discrete foliation. Each leaf - labeled by a discrete lattice time t n - is a piecewise 
linear spatial hypersurface of a given, fixed topology T . More precisely, each 
spatial hypersurface is a simplicial manifold built out of equilateral flat tetrahedra 
with link length a s , together forming a three-dimensional triangulation. Adjacent 
spatial hypersurfaces are then connected with the help of four-simplices, in such a 
way that the whole spacetime becomes a four-dimensional simplicial manifold of 
topology T x [0, 1]. In practice we use T = S 3 and furthermore impose periodic 
boundary conditions in the time direction, such that the spacetime topology 
becomes S 3 x S 1 , but this should be regarded merely as a choice of convenience. 

The resulting four-dimensional triangulation forms a simplicial complex, con- 
sisting of vertices, links, triangles, tetrahedra and four-simplices. The Lorentzian 
nature of the metric properties of the triangulations implies a more refined cat- 
egorization of their elementary building blocks. There are two different types of 
links, depending on whether their endpoints lie on the same hypersurface (lead- 
ing to spatial links) or on two neighbouring hypersurfaces (resulting in timelike 
links). In CDT quantum gravity we allow these two types to have different abso- 
lute length and define the (squared) length of the timelike links to be a\ = aa 2 s , 
with a relative scaling parameter a < 0. Consequently, higher- dimensional (sub-) 
simplices of a CDT geometry also come in different varieties, depending on the 
space- or timelike character of their one-dimensional subsimplices. 

The path integral now becomes a sum over all geometrically inequivalent 
triangulations of this type with a fixed number of time steps. As an action 
we use the Einstein-Hilbert action, which has a natural realization on piecewise 
linear geometries, the so-called Regge action [20J. The corresponding model in 
two dimensions can be solved analytically [21J, but already the three-dimensional 
model has only been solved partially and for restricted classes of triangulations 
[2"2"j 123] • In four dimensions analytical methods are mostly unavailable and one 
must resort to Monte Carlo simulations to extract physical results. To do so one 
needs to convert the path integral into a statistical partition function by applying 
a Wick rotation. 

Because of the presence of a foliation and an associated global notion of time 
one can perform a Wick rotation at the level of the geometries by simply rotating 
a —7- —a in the lower-half complex plane. The Regge action changes accordingly 
and becomes the Euclidean Regge action 



where N 0l N4 and AQ ' denote the numbers of vertices, four-simplices and four- 
simplices of type (4,1) (with four vertices on one hypersurface and the fifth on a 
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Figure 1: Visualization of the phase diagram of CDT quantum gravity, based on 
actual measurements, illustrating the distinct volume profiles characterizing the 
three observed phases A, B and C. 

neighbouring hypersurface). The three couplings k,q, ^4 and A can be expressed 
as functions of the bare gravitational coupling, the bare cosmological coupling 
and the asymmetry parameter a introduced above. By defining kI = K4 + A one 
obtains the version of the Euclidean Regge action which is ultimately implemented 
in the computer simulations, namely, 

S Regge = -koNo + k±N 4 + A(Ni iA) - 6N ), (2) 

where we have grouped together the terms proportional to A. For later conve- 
nience we introduce the notation conj (A) :=N^ 1] -6N for the quantity conjugate 
to A in the action. 

In this way the path integral of quantum gravity is turned into a statistical 
partition function. We will use the freedom to switch to a different ensemble, 
obtained by holding fixed N4 - measuring the size (the total discrete spacetime 
volume) of the system - instead of its conjugate K4. We therefore treat N4 as 
a finite-size scaling parameter which does not appear in the phase diagram of 
the putative continuum theory. The remaining two couplings k,q and A span the 
phase diagram of the system, which we will explore in the next section. 
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3 The phase diagram of CDT 



The first, qualitative description of the CDT phase diagram appeared in [TT] . 
followed more recently by a quantitative and more detailed analysis in [23]. The 
new visualization of the phase diagram shown in Fig. [T]is based on phase transition 
data from the latter. The phase diagram contains three phases, labelled A, B and 
C [T7]. They can be distinguished by looking at a particular feature of the large- 
scale geometry of their ground state, the spatial volume profile N^(t), which 
measures the three-volume in lattice units as a function of proper time t. 

Representative profiles have been visualized in Fig. H]by rotating the function 
N 3 (t) around the t-axis, thereby creating a body of revolution. The average geo- 
metry ("average" in the sense of expectation values) found in phase C exhibits 
the scaling behaviour of a genuinely four- dimensional universe, whose average 
volume profile is consistent with a Euclidean de Sitter spacetime [T5l [T9] . The 
situation in the other phases is very different. The typical volume profile of 
a configuration in phase A consists of an essentially uncorrelated sequence of 
spatial slices, while the configurations in phase B are characterized by an almost 
vanishing time extension, in the sense that almost the entire volume of the system 
is concentrated around a single spatial slice. 

We should mention at this stage that several qualitative features of the CDT 
phase diagram bear a striking resemblance to those of a similar diagram in so- 
called Hofava-Lifshitz gravity [25J, as has been pointed out elsewhere [251 [23]- 
Trying to understand whether this indicates a deeper relation between the two 
quantum-gravitational frameworks provides an additional incentive for perform- 
ing a more detailed quantitative analysis of the phase diagram, as we are doing 
here. 

Fig. [2] shows the quantitative phase diagram, measured at system size N± = 
80.000. The location of the phase transitions depends on the system size: when N4 
is increased, the A-C transition moves towards larger values of kq, while the B-C 
transition moves towards larger values of A. Using finite-size scaling techniques it 
is possible to determine the location of transition points for infinite four-volume. 

Before embarking on a detailed analysis of the CDT phase transitions, let us 
comment on some technical issues with regard to measuring the location of a phase 
transition. The very notion of a phase transition point is of course ambiguous for 
finite systems, because strictly speaking phase transitions can only occur in the 
limit of infinite size. For example, consider the susceptibility xo — (O 2 ) ~ (C) 2 
associated with some observable O. Plotting the function Xo{ K o — const, A) near 
the B-C transition, one will find that it has a maximum at some value A c . The 
location of this maximum can be used as a definition of the transition point, but 
one must keep in mind that its precise value will in general depend weakly on O. 

Another definition of phase transition comes from studying histograms of 
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Figure 2: Phase diagram of CDT quantum gravity [21]. Crosses represent actual 
measurements, dashed lines are extrapolations. 

suitable observables. Consider the Monte Carlo time evolution of the quantity 
No near the A-C transition, as depicted in Fig. [3] (top left). We observe that iVo 
for some time fluctuates around one value, then makes a transition to another 
value, around which it fluctuates for a while, before flipping back to the first value 
again. This behaviour is characteristic of a phase transition, with the two different 
states corresponding to the two phases involved. The associated histogram of Nq 
(Fig. [31 bottom left) shows a double-peak structure. Varying k while holding 
A fixed will change the relative heights of the two peaks. We can now fine-tune 

so that both peaks have the same height. This gives us another definition of 
the phase transition point, although often not the most practical one, since the 
fine-tuning can be very time-consuming. Also, there may be regions of coupling 
constant space or volume sizes where the flipping between the two states is not 
strong enough, in which case the histogram will not exhibit two peaks, but look 
more like a deformed Gaussian distribution. 

As we move towards N^ — > oo, all of these definitions will converge to one 
and the same phase transition point at infinite four-volume. However, what is 
important for our purposes is the observation that already for moderately large 
system sizes those definitions will produce values so close to each other as to be 
practically indistinguishable on a plot like that of Fig. [21 In order to determine the 
phase diagram it is therefore sufficient to fine-tune the parameters until a flipping 
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Figure 3: Monte Carlo time evolution (i) of iVo at the A-C transition (top left) 
with associated histogram (bottom left), (ii) of conj(A) at the B-C transition 
(top right) with associated histogram (bottom right). 



of phases is observed. Not all observables are equally suited for this kind of 
analysis, because the flipping between phases can occur with different amplitudes 
for different observables. A good choice is usually the quantity conjugate to the 
coupling that needs to be fine-tuned. The form of the action (j2J) suggests using 
the quantity No for the A-C transition, and conj (A) = N^ 1} -6N for the B-C 
transition. The Monte Carlo time evolution and associated histogram for the 
latter are shown in Fig. [3] on the right (top and bottom). 

It is instructive to track how the phase transition signature changes as one 
moves along the phase transition lines, while holding the system size iV 4 = 80.000 
fixed. Along the A-C transition we have not observed any appreciable change, 
although we have not investigated the immediate neighbourhood of the triple 
point where all three phases meet. The situation for the B-C transition is rather 
different. Let us start on the B-C line at the point k,q = 2.2 and move in both 
directions along the line. As we move to the left, the jump in conj (A) associated 
with the phase flip decreases, and around k = 1.0 no observable signature of 
a phase transition remains. We conclude that the B-C transition line has an 
endpoint, which for N± = 80.000 is located around k = 1.0. Its precise location 
remains to be determined. 

Conversely, when moving to the right from k = 2.2, we observe that the 
jump in conj (A) grows rather quickly, until at kq = 2.5 we no longer observe the 
phase flipping, although the reason for this is very different from the situation 
encountered at the left endpoint. The increased size of the jump in conj (A) is 
accompanied by a deepening of the vertical gap separating the two peaks of the 
histogram of conj (A). As it becomes deeper, the frequency of the phase flipping 
decreases significantly, until it becomes so small that no flipping is observed 
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during the entire simulation. This indicates that the simulation becomes stuck 
in a metastable state and does not reach thermal equilibrium. In Fig. [2] this 
region is marked by a dashed line; here conventional methods are insufficient to 
measure the phase transition location with acceptable accuracy. We are currently 
investigating the use of multicanonical Monte Carlo simulations [27], which at 
least in principle can offer a solution to this problem. 

4 How to measure the order of phase transitions 

Measuring the order of phase transitions requires some care, as is illustrated 
by the history of dynamical triangulations. The quantum gravity model based 
on four-dimensional Euclidean triangulations (which lack the foliation present in 
CDT) has a phase transition, which initially was determined to be of second order 
[28l [29j [30] . However, later studies with larger systems established instead that 
the phase transition is of first order [311 [32] . 

Criteria to distinguish between first- and second-order transitions have been 
summarized conveniently in [33], although care should be taken in extrapolating 
from conventional systems, simulated on static lattices, to dynamical triangu- 
lations, where the lattice itself is dynamical. Let us begin by considering the 
situation of infinite system size. In this case first-order transitions are uniquely 
characterized by the existence of observables (first derivatives of the free energy) 
that are discontinuous across the transition. In addition, the size of the fluctua- 
tions relative to the average tends to zero in the infinite- volume limit, implying 
that the probability distribution of such an observable approaches the sum of two 
delta-function distributions. The displacement between the two peaks is precisely 
the difference between the two equilibrium values of the observable. 

The probability distribution of an observable in the infinite- volume limit can 
be understood as arising through a limiting process of the corresponding distribu- 
tions for finite volume. Their measured counterparts are histograms of the kind 
displayed in Fig. [31 It is sometimes stated in the literature that a double-peak 
structure of a histogram signals a first-order transition, but this is somewhat mis- 
leading and can even be wrong. Rather, the existence of a double-peak histogram 
allows one to confirm the first-order nature of a transition by considering a se- 
quence of histograms for increasing system size. We are dealing with a first-order 
transition whenever the double-peak structure becomes more pronounced with 
increasing volume. Quantitatively, this means that the vertical gap associated 
with the double peak increases. By "vertical gap" we mean the difference be- 
tween the peak heights (assuming they are equal) and the height of the minimum 
in between them. How the vertical gap changes with system size is the method of 
choice to confirm the first-order nature of a transition, provided one can simulate 
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sufficiently large systems, which develop a double-peak structure. One could in 
principle try to use the same method to conclude that a transition is of second 
or higher order by measuring the vertical gap as a function of inverse system size 
and extrapolating to the limit of infinite volume, where it should vanish. How- 
ever, because it is often difficult to measure this gap with good accuracy, showing 
conclusively that it really goes to zero in the limit can be problematic. Fortu- 
nately, there are better criteria at hand - involving the measurement of critical 
exponents - which one can use to establish that a given transition is of higher 
order. 

One such exponent measures the shift of a transition point with system size. 
Recall first how this works for a conventional lattice system such as the Ising 
model with volume V = L d , where d is the system's dimension [31]. Considering 
the temperature-driven phase transition of the Ising model and using the location 
of the maximum of the magnetic susceptibility to define a transition point (3 C (V), 
one finds a power-law behaviour 



for sufficiently large system size. The exponent v governs the increase of the 
correlation length in a second-order transition as one moves towards the critical 
point /3 c (oo) on an infinite lattice. For first-order transitions there is no correlation 
length and one expects the specific scaling [33] 



A sufficiently strong violation of v = 1 therefore signals the presence of a second- 
order transition. 

To work with the criteria (J3J or (jl]) one needs some way of judging whether the 
system sizes under consideration are large enough. Given N data points ordered 
by system size, we can make M fits where the i-th fit, i G {1, . . . , M}, is made by 
using the restricted set of data points with labels i, . . . ,N. If the corresponding 
sequence of exponents is drifting, it can be an indication that the system sizes are 
not large enough. Obviously, this method is only useful when one has sufficiently 
many data points. 

Another quantity of interest is the so-called Binder cumulant associated with 
an observable O, which may be defined as [33J 



and is always nonpositive. Considering B as a function of the couplings, its 
local minima are at transition points, where the fluctuations are largest. We can 
measure these minima for different system sizes and by extrapolation determine 




(3) 




(4) 
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Bq 111 (1/N a = 0). This quantity is zero when the probability distribution of O 
approaches a delta function around an expectation value (O) in the infinite- 
volume limit, as is expected at a second-order transition. On the other hand, 
many aspects of first-order transitions are described well by approximating the 
histogram of the observable with a superposition of two distributions centred at 
the expectation values (0) — 1 and (0} = 2 in the two phases. Again, if these 
distributions approach delta functions in the infinite- volume limit, the minimum 
of Bo will be given by 

(O 2 - O 2 ) 2 

^(1/^ = 0) = -^^, (6) 

which is obtained for a value of the coupling constants where the relative strength 
of the delta functions at d and 2 is given by 2 2 /{0\ + 2 2 ) and 0\J{0\ + Of) 
respectively. We conclude that a sufficiently strong deviation of Bq 111 from zero 
signals a first-order phase transition. Using the value of the Binder cumulant 
to make the case for a second-order transition is more difficult, because a weak 
first-order transition may show a convergence to a value close to zero. 

There are other methods, involving the measurement of other critical expo- 
nents, which can help to determine the order of a phase transition, but their 
translation to CDT quantum gravity leads to ambiguities. We will therefore pro- 
ceed by applying the methods described above to the A-C and B-C transitions 
found in the CDT system. The third transition (A-B) is currently of minor inter- 
est, since it bounds two phases that most likely have no relevance for continuum 
physics. 



4.1 The order of the A-C transition 

The simulations for the analysis of the A-C transition were run at A = 0.6, with 
system sizes ranging from 40/c to 150/c. The data presented in the following are 
the results of simulations where the number of sweeps was approximately 5 • 10 6 , 
with one sweep representing one million attempted Monte Carlo moves. 

Our simulations suffer from a slow convergence of observables, which means 
that the standard error algorithms produce error values which significantly under- 
estimate the true uncertainties of the measurements. In all our measurements we 
have attempted to obtain more realistic error values by systematically studying 
the convergence of observables. Given a set of data samples at Monte Carlo times 
Ti, . . . , r n , we can define a time- dependent average (O) (f ) by including only those 
data samples for which T\ < r < f . We then plot this quantity as a function of 
f and try to extract reasonable error values. This method introduces a degree of 
subjectivity, but in our opinion produces more realistic error estimates. 
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Figure 4: Histograms of N at the A-C transition for four different system sizes. 
The histograms have been normalized to obtain probability distributions P. 



We first perform a histogram analysis. Fig. H] shows the histograms for the ob- 
servable N Q (more precisely, the rescaled quantity Nq/N^) at four different system 
sizes. The histogram at = A0k assumes an approximate Gaussian shape, which 
gets distorted at N± = 80k, with a double peak starting to emerge at N4 = 100k, 
which becomes more pronounced at = 120k. For the two larger volumes, the 
vertical gap associated with the double-peak clearly increases, while the mutual 
distance of the peaks stays approximately the same. This behaviour is a clear 
signature of a first-order transition. The emergence of the double-peak shape 
seems to occur somewhere between N 4 = 80k and N4 = 100k. This may be com- 
pared with the analysis performed in the context of four-dimensional Euclidean 
dynamical triangulations, where the double peak emerged between N4 = lQk and 
N4 = 32k [31]. The more rigid structure imposed on the triangulations by the 
causality conditions of CDT may be responsible for the shift of the appearance 
of the first-order signal to somewhat larger volumes. 

Let us try to find additional evidence for the first-order nature of the A- 
C transition by first measuring the shift exponent v defined in eq. (jl]). For 
sufficiently large system sizes we expect a power-law behaviour 

k c (N a ) = ^(oo) - CN^'\ (7) 

where KqIN^) denotes the transition point at system size N4 and C is a propor- 
tionality factor. 
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Figure 5: Measured A-C transition points Kq(N^) at A = 0.6 for different system 
sizes A/4, together with a best fit to the inverse- volume expansion eq. (jHJ). 

We have already mentioned that the notion of a transition point is ambiguous 
for finite system sizes. For our present purposes, we define Kq(N±) as the location 
of the maximum of the susceptibility xa^o = C^o) — (^o) 2 - To measure this 
maximum we use an extrapolation method due to Ferrenberg and Swendsen |35j, 
also known as histogram method [51] . 

The best fit using all data points gives v =1.11(2), with a reduced chi-squared 
of x 2 = 2.1. Unfortunately, the three-parameter fit ([7j) is not very stable when 
removing the points with the lowest values of N4. Taking this into account, it 
becomes clear that the error bar of 0.02 significantly underestimates the error 
associated with the determination of v from the data. It indicates the presence 
of subleading terms which for the volume range studied interfere with the pure 
power-law behaviour of (|7j). Consequently, the strongest conclusion we can draw 
from using this method is that the measured z> = 1.11 is compatible with the value 
z> = 1 characteristic of a first-order transition. 

As an additional cross-check we have made another three-parameter fit to the 
generic functional form of a large- A/4 expansion of Kq(N^), valid at a first-order 
transition [53], namely, 

k c (N 4 ) = ^(oo) - C/N 4 - D/Nl + 0(l/Nl). (8) 
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Figure 6: Dependence of the minimum -B™™ of the Binder cumulant B^ on the 
system size iV 4 , at the A-C transition and for A = 0.6. 



Fig.[5]shows the measured data points and the best fit, based on the expansion (JH]). 
The errors associated with individual data points are too small to be visualized in 
a sensible way. The quality of this fit is comparable with that obtained from fitting 
to the power law ([7]). The contribution from the l/iVf-term is subdominant to 
that from the l/A^-term, as should be. To summarize, we can say with confidence 
that the shift data are fully compatible with a first-order transition. However, it 
is also clear that if one wanted to go beyond this statement and nail down the 
value of the critical exponent v with better precision, one would need to consider 
larger system sizes. 

Our next and final step in analyzing the A-C transition will be to consider 
the Binder cumulant for N , 



B 



N 




(9) 



N at the A-C phase transition. Like in the case of 



which has a local minimum B 
the shift exponent u, we use the histogram method to determine this minimum. 
In Fig. |6] we have plotted B™™ as a function of system size. The error bars are 
much larger than for the measurement of the shift exponent, where we measured 
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Figure 7: Histograms of conj(A) at the B-C transition for three different system 
sizes. The histograms have been normalized to obtain probability distributions 
P. 

the location of the susceptibility maximum, and not the maximum value itself. 
The plot shows clearly that moves away from zero as one goes to large 

system sizes, reconfirming the first-order nature of the transition. 

By contrast, when going to smaller volumes one finds a different behaviour, 
where -B™™ increases before reaching a local maximum. Evidently, for these small 
volumes the system lies outside the scaling region, where quantities like B™" 1 are 
expected to behave according to a power law. Looking at the plot, a rough 
estimate of the onset of the scaling region seems to be around or above A4 = 60k. 
Comparing with Fig. HJ this may be correlated with the observed emergence of 
the double-peak structure in the histograms. 

4.2 The order of the B-C transition 

Continuing our investigation of the order of phase transitions in CDT quantum 
gravity, we now turn to the B-C transition. In this case, we have fixed the inverse 
gravitational coupling to k = 2.2 and analysed the system at sizes between 40k 
and lQOk. The number of sweeps used was approximately 2.5 -10 6 , with one sweep 
again corresponding to one million attempted Monte Carlo moves. Fig. [7] shows 
the histograms of the quantity conj(A) = A^ 4 ' 1 ^ — QN for three different system 
sizes. The situation differs from that of the A-C transition in that we observe a 
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Figure 8: Measured B-C transition points A c (iV4) at k, q = 2.2 for different system 
sizes iV 4 , with a best fit to eq. ([TO]) to determine the shift exponent v. 

double-peak structure for all four-volumes. The peaks are most pronounced for 
smaller volume and appear to be merging when the volume is increased. The 
plot shows that their mutual distance decreases, roughly like I/-/V4. This is a first 
indication that the B-C transition may not be a first-order transition^] 

We proceed by measuring the shift exponent for A. Analogous to what we 
did for the A-C transition, we will use the formula 

A c (iV 4 ) = A c (oo) - CN^ 1/0 , (10) 

where A c is defined as the location of the maximum of the susceptibility Xconj(A) = 
(conj(A) 2 )— (conj(A)) 2 . Data points and best fit are displayed in Fig. [SJ We have 
not included any error bars because they turned out to be too small. The best 
fit through all data points yields v = 2.39(3). To judge whether our range of 
system sizes lies inside the scaling region we have again performed a sequence of 
fits by successively removing the data points with the lowest four-volume. The 
corresponding values for v are 2.39(3), 2.51(3), 2.49(3) and 2.51(5), where the 

1 Since we have not performed a time-consuming fine-tuning of A to obtain peaks of equal 
height, we have not been able to extract a reliable estimate of their vertical gap (relative height 
of peaks w.r.t. minimum in between peaks, as defined earlier in Sec.|l]). 
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Figure 9: Dependence of the minimum -B^nj(A) °f the Binder cumulant -B CO nj(A) 
on the (inverse) system size , at the B-C transition and for kq = 2.2. (Fit 
excludes the two points on the right.) 

last fit was done with all but five data points removed. Again the error bars are 
based on making cuts in the sampled data as described above when discussing 
the A-C transition. The sequence suggests that the data point with the lowest 
four-volume lies outside the scaling region. Removing it from the fit we get 

P = 2.51(3). (11) 

This result makes a strong case for a second-order transition, since the prediction 
v = 1 for a first-order transition is clearly violated. 

Let us finally consider how the minimum of the Binder cumulant Bq of def- 
inition (jSJ), with O = conj(A), depends on the system size. In Fig. ED we show 
(A) as a function of inverse system size. The errors turned out to be much 
smaller than for the corresponding measurements at the A-C transition, with er- 
ror magnitude approximately equal to the radius of the dots in the plot. Inside 
the scaling region the minimum of the Binder cumulant is expected to behave like 
a power-law function. To determine whether our data points lie inside this scaling 
region, we have again performed a sequence of fits by successively removing the 
data points with the lowest four-volume. We have found that the data points at 
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Observable O 


B$ a (N 4 -> oo) 


conj(A) 


-0.003(4) 




-0.001(3) 


N 2 


-0.000 0001(3) 


Ni 


-0.000 003(7) 


N 


0.000 0(3) 



Table 1: Measurements of 5™ n (iV4 — > oo) for various observables O, where Nk 
denotes the number of fc-dimensional (sub-)simplices in the triangulation. 

N 4 = 40k and N 4 = 50/c both lie outside the scaling region. The curve displayed 
in Fig. corresponds to the fit with those two points removed. 

Table [1] collects the results of measuring BQ™{N i — > oo) for five different 
observables O. All of these observables exhibit a phase flipping at the transition, 
with associated double-peak histograms. The interpolation to N 4 — > oo is made 
by assuming that the quantities B^ in {N 4 ) for sufficiently large volumes have a 
scaling behaviour like in eqs. ((7j) and (flOl) . namely, 

B™ n {N 4 ) = B™ n (N A oo) - C{0)/Nl {0 \ (12) 

with (9-dependent exponents p{0) and coefficients C(0). If we had to rely on 
the Binder cumulant measurements alone to determine the order of the B-C tran- 
sition, we would encounter the already mentioned problem of how to distinguish 
between a second-order transition and a first-order transition where Bq 111 con- 
verges to a value close to zero. In the case at hand, in addition to the evidence 
already presented in favour of a second-order transition, the size and quality of 
the measurement errors for the minimum of the Binder cumulants are such that 
we can conclude with confidence that the results displayed in Table [Hare certainly 
consistent with a limiting value of zero for ^4—7-00. This would be true even if 
the errors overestimated the true uncertainties by a factor of two, say. 

5 Discussion and outlook 

We set out with the goal to determine the order of the two physically relevant 
phase transitions in CDT quantum gravity. For conventional systems on static 
lattices, describing physics on a fixed background geometry, the methods for doing 
so are well established. It is not clear a priori whether all of them are applicable 
for systems based on dynamical lattices, such as CDT, which reflect the dynam- 
ical character of the quantum-gravitational, geometric degrees of freedom they 
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aim to describe. Evidence that standard methods can be adapted rather straight- 
forwardly to systems of dynamical geometry comes from Euclidean dynamically 
triangulated (DT) systems. This has been demonstrated for two-dimensional 
Euclidean quantum gravity (coupled to matter), where extensive computer simu- 
lations [36] are in agreement with Liouville quantum gravity, which can be solved 
analytically. It is also true in four-dimensional DT. Although this model does 
not appear to describe four-dimensional quantum gravity, it does exhibit a well- 
defined scaling behaviour in the so-called branched-polymer phase and finite-size 
scaling works very well [30J. 




However, since causal dynamical triangulations differ from their purely Eu- 
clidean counterparts both in their set-up and (hopefully) the continuum physics 
they describe, we decided to not take the validity of these methods for granted. 
Rather, we used several independent criteria for determining the transition orders, 
and checked whether the results were mutually consistent. 

Our analysis consisted of three parts: the histogram analysis, the measure- 
ment of the shift exponent and the analysis of Binder cumulants. The histogram 
analysis was sufficient to confirm the first-order nature of the A-C transition. It 
probes the definition of a first-order transition in a direct way and is therefore the 
method of choice to confirm the first-order nature of a transition. Measurements 
of the Binder cumulants also pointed rather unambiguously to a first-order transi- 
tion. The absolute values of the minima were increasing for increasing volume N4, 
contrary to what one would expect for a second-order transition. Measurements 
of the shift exponent, which governs the dependence of transition points on the 
system size, were less conclusive. The results were compatible with v — 1, valid 
at a first-order transition. However, since the stability of the curve fits under 
successive removal of data points with small volume was not particularly good, 
we did not succeed in determining the exponent v convincingly from the data, 
but only in establishing its consistency with the first-order value. 

The corresponding analysis of the B-C transition gave a very different picture. 
The histogram analysis showed a double-peak structure, but without the charac- 
teristic features of a first-order transition. The distance between the two peaks 
diminished with increasing N4, indicating a higher-order transition. This was 
corroborated by measuring the shift exponent, whose value v = 2.51(3) entails a 
strong violation of the prediction v — 1 for a first-order transition. In addition, 
the results of the Binder cumulant analysis were also clearly consistent with the 
presence of a second-order B-C transition, as announced previously in [16] . 

In summary, we have found strong evidence that in CDT quantum gravity 
the A-C transition line is first order, while the B-C transition line is second 
order. The latter result is both remarkable and attractive. It opens the door 
to studying critical phenomena in CDT and to defining a continuum limit for 
vanishing lattice spacing (UV cutoff). A discussion of the effective gravitational 
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coupling constant, and its dependence on the bare couplings of the model - inside 
phase C - was initiated in [T5]. In view of the results presented here, a future task 
will be to extend this analysis by studying the flow of this coupling constant when 
approaching the second-order B-C phase transition line. Apart from the intrinsic 
interest in understanding the behaviour of quantum-gravitational observables in 
this limit, which appears to be associated with probing (sub-)Planckian physics, 
this may enable one to make a more explicit connection with the asymptotic 
safety scenario, as outlined recently in [13], as well as anisotropic gravity models 
of Hofava-Lifshitz type. 

To make this connection more concrete, let us consider a potential application 
of our analysis, based on the presence of a second-order transition. As shown in 
[T4"l[T7] . well inside phase C the time extent T of the universe scales like iV 4 1//4 , and 
the three- volume V^it) of a generic spatial slice at time t like N^ 4 . By contrast, 
in phase B this time extent essentially vanishes, since all three-volume is located 
at a single spatial slice. Approaching the B-C transition line from inside phase 
C, we observe that T, measured in lattice time steps, decreases, at least in the 
region k < 2.4 where detailed measurements are available. Its precise behaviour 
still needs to be determined, but it is natural to conjecture a behaviour of the 
form 

T C (N 4 ) oc (A c (oo) - A c (iV 4 )) e iV 4 1/4 (13) 

at the pseudo-critical point A c (A r 4). Combining this with the scaling ( flUj) of the 
coupling A, we obtain 

T c (iV 4 ) = Nl /4 ~ e/i> , (14) 

implying that the time extent T scales anomalously. Interestingly, such an anoma- 
lous scaling of time relative to space in the ultraviolet is also a key feature of 
Hofava-Lifshitz gravities [25]. Following this line of reasoning further, one could 
imagine that both v and 9 depended on kq, in this way potentially allowing for 
both anisotropic and isotropic UV completions, depending on where the B-C 
phase transition line is approached. 
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